Tricritical behaviour of Ising spin glasses with charge fluctuations 
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We show that tricritical points displaying unusal behaviour exist in phase diagrams of fermionic 
Ising spin glasses as the chemical potential or the filling assume characteristic values. Exact results 
for infinite range interaction and a one loop renormalization group analysis of thermal tricritical 
fluctuations for finite range models are presented. Surprising similarities with zero temperature 
transitions and a new T = tricritical point of metallic quantum spin glasses are derived. 
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Spin glass phase transitions can have important effects 
on a variety of characteristics of fermionic systems at low 
temperatures including nonmagnetic properties, and vice 
versa. This has been emphasized in a series of recent the- 
oretical articles on quantum spin glass transitions [Q-^. 
Considerable interest in these problems was raised by ex- 
perimental results for heavy fermion systems |^J^. Spin 
glass phases known to exist between antiferromagnetic- 
and superconducting phases of LaSrCuO phase diagrams 
1^ may also be considered under the aspect of interplay 
of magnetic and electronic properties. 

In this Letter we report results on existence and prop- 
erties of tricritical fermionic spin glass transitions which 
happen to emerge as the fermion concentration, rele- 
vant for the effective spin dilution, is moved through a 
characteristic value. We analyzed in detail the tricriti- 
cal behaviour of the Ising spin glass (denoted by ISGj) 
on fermionic space with 4 states per site instead of the 
usual two of SK models for example. Results are given 
i) for infinite-range-, and (by use of the renormalization 
group) ii) for finite-range spin interaction, and iii) for a 
metallic model with additional electron hopping hamilto- 
nian. The two nonmagnetic states per site provided by 
the fermionic space allow the system to adjust its effec- 
tive random spin dilution according to quantum statis- 
tics. On decreasing the effective spin density and hence 
Tc the system is driven through a tricritical point into a 
regime of discontinuous spin glass transitions. The tri- 
critical point {TCP) turns out to be particularly inter- 
esting, since quantities such as density of states, fermion 
concentration (for given /i and vice versa), local suscep- 
tibility, and spin correlation function behave nonanalyt- 
ically at the tricritical spin glass transition, and thereby 
change substantially critical properties of those quanti- 
ties which typically define and signal spin glass transi- 
tions. Experimental observation is also favoured by sev- 
eral aspects such as increased upper critical dimension 
di"'* = 8 (instead of a decreasing from 4 to 3 in 
conventional 0^-thcories [^), which allows critical fluc- 
tuations in 31? to have much stronger effects than the 
usual logarithmic corrections and the occurence of phase 



separation in the 1st order regime, a phenomenon shared 
by the phase diagram of the BEG-modcl for He^ — He'^ 
mixture s pj ] . We also compare with other classical spin-1 
models 0. 

Understanding the key features of the ISGf as the 
simplest model which takes spin glass order and charge 
correlations into account provides a useful guide to phase 
diagrams of spin glasses allowing for thermally activated 
hopping or metallic conductivity. It appears to be generic 
for the behaviour of an even larger class of models, in a 
way comparable with the BEG-model The ISGf is 
defined by the hamiltonian H = —-^^ Jijcricrj — iJ^Ui 
with spins a = n| — n^, the fermion number operator 
n = n| + rt| , and gaussian-distributed with vari- 
ance J. All spin- or charge-correlations of this model 
are static. As an example the local susceptibility Xai'^) 
has the Fourier transform x(ci') — I3{q~ q{x))6i^o with 
q —< (Ji[t)ai{t') > and q{x) denoting the Parisi solution 
of the ISG f . In contrast to the standard 2-states per site 
Ising spin glass this 4-states per site model feels quantum 
statistics due to the relative occupation of magnetic and 
nonmagnetic states. Quantum dynamics however enters 
in the electron Greens function and correlations being 
defined with an odd number of equal-time fermion op- 
erators. The fermionic path integral representation not 
only produces the correct spin-static field theory for the 
spin glass transition, it also helped us to find a close rela- 
tionship between nonanalytical behaviour of the density 
of states or fermion concentration and the special features 
of the tricritical spin glass transition. In the same way 
we identify the surprising similarity between this classical 
tricritical theory for finite Tc and the recently analyzed 
quantum theory of the T = 0-quantum paramagnet to 
spin glass transition of a metallic model. 

Working in the grand canonical ensemble our de- 
tailed analysis of the ISGf shows that upon increas- 
ing the chemical potential from /j, — at half-filling 
the line of continuous spin-glass transitions given by 
Tc — Jq{Tc) is only realized up to a tricritical point lo- 
cated at Tc3 ~ J/3. Beyond this point, ie for Tc < ^ 
one enters the domain of discontinuous spin-glass tran- 
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sitions. Fig. ^ displays this most interesting part of the 
phase diagram which surrounds the TCP. All relevant 
zero-field properties can be concluded from our result 
for the saddle-point free energy and expanded around 
the TCP, which reads (the reader may also extract the 
tricritical behaviour from the simpler approximate form 
using q[x) — qea constant) 

3 3 

/ - frcP = - Aic3 + Jli^^gg - rrT^)6q + -T^rliSqf 

dxq^{x)] - Tr[ / dx[xq^{x) + "iqix) I dyq^{y)] 

^ Jo JQ 

3Sq£ dxq^x) + ^{Sqf] - ^ dxq\x)}, (1) 



. The constants 

3 A-2\ 



where Sq = q - qTCP,gJ = A* - Mcs 
as nonordering field, and t = '^Z,'^"^ 

^ c3 

are given by rg — ^,rr = 2(1 — |C with ( 
tanh{^c3,/Tc'i) — 0.9938, and /ic3 — ^arcosh{2exp{^)) ~ 
0.9611J as the characteristic chemical potential locating 
the TCP. The average filling factor corresponding to /ic3 
is evaluated as < >— 1.6625. By symmetry w.r.t. 
fj, — one also finds a TCP for less than half-filling at 
< j/c3 >— 0.3375. Thus the low- and the high-filling do- 
mains host discontinuous spin glass transitions. We shall 
argue below that this is a rather general phenomenon. 




FIG. 1. Vicinity of tricritical point (TCP) for positive 
chemical potential. Continuous spin glass transitions occur 
on curve (a) above the TCP (thick unbroken line). Below 
the TCP first order thermodynamic transitions take place on 
curve (c). Curve (d), starting at the TCP, and curve (b) 
below the TCP limit the existence regime of ordered and dis- 
ordered phases respectively. 

The particularity of the TCP stems from the fact that 
the replica-diagonal fields Q"" become critical in addi- 
tion to the usual replica-overlap fields Q"** with a ^ b. 
On approaching the TCP with T — s- Tcs and fi fixed for 
example, the saddle-point solution q = Qgp, as obtained 
from Eq. dl]) , develops nonanalytical behaviour given by 



Sq = q{T) - qiT,^) = A±{±t)^ + 0(r) 



(2) 



with A 



\/3A_ ~ 0.1999. This behaviour is at 
the origin of the crossover from the critical exponents 



—a = f3q = ^ = \, vahd for Tc > Tc3, to the tricritical 
ones as = /3g3 = 73 = ^ obtained here for T Tcs. 
The free energy given above can be cast into the scal- 
ing form / = t'^~"- G{^) for the disordered phase with 
a — —\ and crossover exponent (j) = 2. In the tricrit- 
ical regime the singular part behaves like f^^^g" ~ 5^ 
whence — from above. In contrast to crystal-field 
split spin glasses [|l^ the quartic coefficient 7/4 of our free 
energy, Eq.(P, is nonzero and one obtains the Parisi so- 
lution q[x) ~ 9^7^ < a; < xi and q{x) = q{l) 
for xi < a; < 1. The plateau height is found to satisfy 
q{l) ~ 5q + 0{5(f). Consequently, plateau and break- 
point scale like \/\t\ + 0{t) at the TCP, while linear 
T-dependence is reserved to T^ > Tcs. Adapting the no- 
tation of we express our result for the irreversible 
response q{l) — q{x) ~ \t\^^ in terms of the exponent 
/3a3 = 1 for r ^ r,3 and /?a = 2 for T ^ > T.g. 
For the Almeida-Thouless line at tricriticality we find 
- ^ tftanh{^))f/^T^J^ + 0{tI^) with 

- — -. Hence we obtain the critical expo- 
while 6* = I for all Tc > T^s. These 
values do not satisfy the scaling relation 6*3 = with 

Pas = 1 + (73 — Q!3)/2. Along the lines described by 
D.Fisher and H.Sompolinsky |ll|] , this problem of mean- 
field exponents will be resolved below by the renormaliza- 
tion group analysis of the coupling j/4 of the finite-range 
and finite-dimensional ISGf. 

The fermionic nature of the ISG / and of related mod- 
els on Fock space calls for a representation in terms 
of electron Green's functions. We shall see that even 
if charge-fields do not occur or do not become critical 
together with the spin-fields, nonanalytic behaviour at 
the spin-glass transition can become not only observable 
but also intimately related to quantities defined in the 
charge-sector. We mention our result for the density of 
states (DOS), which is given in the disordered phase by 



TT^ 8lV3V-^ 

_ Tc3-Ta 

TAT = 7^ 



nent 63 = ^ near Tc3 



< Pa{e) > = 



ch{l3y) + ch{l3{e + fi)) 
/2^J[ch{f3fi) + exp{^^)] 



(3) 
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For T < Tc we note that < p >=< p > \q=o + 0{q^) 
and the fact that the self-consistency equations can 
be obtained from Eq.(||). The internal energy U — 
'ZJ deef{e) < p„[e) > is found as U = -1.996J - 
1.045Sq + 0[5(f,q'^). Inferring Eq.(||) the specific heat 
is seen to diverge like c ^ confirming thus ^3 = i 

on both sides of the TCP. Transitions at Tc > Tcs 
obey Sq ~ g(l) ~ t instead. The average filling fac- 
tor < ly >= J2cr < n„ > obeys the relation < v 
1 + tanh{f3fi){l — q) (we find that this relation is invari- 
ant under replica symmetry breaking). For fixed p, this 
implies that < > shows nonanalytical -^/[rl-behaviour 
near the TCP, and the same holds for the electronic den- 
sity of states. Collecting the results by 
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d < p^je) > d<v> dq_ I |_i/2 

dT ^ dT ^ dT^^ ^ ^ ^ 

we conclude that all these quantities diverge with the 
MF-exponent 03 of the specific heat both from above 
and below Tcs- The local susceptibility x = JiJ+Pji^ + 

6q — <lix)) has a divergent slope for T | Tcs but a fi- 
nite one for T f Tcs due to our result q{l) = Sq + 0{6q^). 
Inferring Eq.(|) one finds dx/dT - t-^6'(t) + 0(t). It 
is thus the nonanalytic change of the fermion concentra- 
tion, which yields in combination with the ubiquituous 
replica-overlap critical fluctuations the special tricritical 
behaviour as the TCP is approached. What happens if 
one wishes to consider fi as a function of given average fill- 
ing < V >? Again the exact relation given above, which 
holds true if p{< v >) is constrained to be nonrandom, 
shows that now the chemical potential will acquire 
corrections. The other conclusions remain unchanged. 

The fermionic picture also allows to study correla- 
tions of the statistical fluctuations of the density of 
states. The DOS'-cumulant < (5p^(e)5p^^''(e') > with 
5p = p— <p>\s sensitive to the spin glass transi- 
tion. We obtain < 5 p"" {0)5 p^ {0) >= ^^^ptlq^ + 0{q^) 
with A(0, 0;0 < \p\ < \pc3\) decreasing monotonuously 
from the maximum value 0.6592 at half filling to 0.3613 
at the tricritical point. Hence this cumulant, being de- 
fined in the charge sector, decays to zero at the spin- 
glass transition quadratically in r for Tc > Tcs and lin- 
early at the TCP. Beyond this point it is discontinuous 
at the transition. The fO^-cumulant also feels replica- 
symmetry breaking and becomes proportional to q^{x). 
In a nonvanishing magnetic field the DOS'-cumulant be- 
comes nonzero also above T = Tc{H = 0). 

As line (b) of Fig.l merges with line (a) at the TCP, 
the replica-diagonal fields Q"" become critical together 
with the off-diagonal fields, a phenomenon very unusual 
for classical thermal spin glass transitions and so far only 
known to occur in special limits |lj] . In addition Q"^"" ap- 
pears linearly in the Lagrangian and hence plays a spe- 
cial role. These crucial features are surprisingly shared 
by the thermal tricritical theory and the T = 0-quantum 
theory for metallic spin glass- [p[ and for the transverse 
field Ising spin glass transitions Q . These crosslinks with 
T = 0-quantum phase transitions with irrelevant quan- 
tum dynamics (for dynamic exponents z > 2) are best 
appreciated in field theory. We derived the Lagrangian 
for the tricritical and finite range ISG / and a Lagrangian 
of the same structure is obtained for generalized models 
(eg with a transport mechanism) at finite temperature 
by integrating out the dynamical degrees of freedom 

L = \j d^xl^ E + ^ E Q""(-V' + u)Q'^- 

1 1 ' 

+ -Tr'i^/Q'^'f - - E Q'"'^'' - y EW"")' 



- yTr'Q3 _ ^ Q-Q-'-Q^" + | ^(Q'^f)^], (5) 

where 4(^)(o) = (^^)(o) = (^^)(o) = f and u^") = de- 
note the bare coefficients at tricriticality. One fourth or- 
der term relevant for replica symmetry breaking is kept. 
Replicas under Y! Tr' are distinct. The Q'^'^Q^^- 
coupling is renormalization group generated as in the 
metallic quantum spin glass and leads to a common upper 
critical dimension di"'' = 8. A shift of the fields remov- 
ing the redundant ((5"'')^-mass term was performed at 
each step of our one-loop renormalization group anal- 
ysis. Since Q"" and Q^^-fields may possess different 
anomalous dimensions, they are kept separately in the 
Lagrangian. Apart from the striking similarities between 
the tricritical ISCf Lagrangian and the metallic quan- 
tum case there are also some important differences: 

i) Time-integrals are absent, each Q-field may be 
viewed as Qa;=o,i^=o(x). 

ii) Unlike the metallic model the linear term is time- 
independent. 

iii) Three relevant cubic couplings appear instead of 
one in the metallic case. 

iv) The dangerously irrelevant (DI) u J dT{Q'^'^)^ quan- 
tum mechanical interaction of the metallic case |5[| turns 
into a relevant mass term for the ISG / . 

While in case of the metallic T = 0-transition the 
M-interaction was seen to render the nonlinear suscep- 
tibility Xni less divergent than the spin glass suscepti- 
bihty XSG, the ISCf shows Xni ^ XSG due to prop- 
erty iv). This was confirmed by our MP calculation of 
the small field behaviour of the magnetization m which 
yields to OiH') m = f - (f (1 - Bfth{^)))i^ 
for T > TcsiH — 0),p — Pc3 near the tricritical point. 
DUE to XsG ~ ^ near > T^g and Xsg ^ W\ 
the Q-propagator l/(fc'^ + m^) yields the MF correla- 
tion length exponent 1^3 = j a.t the TCP and = 5 
in the continuous regime above. As for the metallic 
T = 0-transition the MF-exponents violate hyperscal- 
ing even in the upper critical dimensions d|,3^ = 8 (TCP 
of ISCf). The ^Q'"'Q''''-term is at the origin of this vi- 
olation for both the tricritical point of the ISCf and 
for the T = 0-metallic spin-glass transition. In the 
metallic case (z — 4) the replacement of dimension d by 
the 'quantum'-mechanical dimension dqm = d + 2z — 9t 
yields modified hyperscaling relations which are satisfied 

(u) 

by MF-exponents in d = dc =8. The same result is 
obtained for the tricritical ISC / by replacing d ^ d— 9t. 
In both cases Ot denotes the dimension of the DIG t. 

We studied tricritical fluctuation in a 1-loop renormal- 
ization group analysis for the Lagrangian (^. Anomalous 
dimensions 77 and 77 are introduced to account for Q"'^''- 

and for Q^'^-fluctuations, respectively. These exponents 

2^2 

and the one for the DIC t are given by rj — ^i^^yi , V — 
2^2 

Trr^,Ot = 2. We obtain the following RG flow equations 
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— = (- - IIkI + 16kiK2 + QKl)r - 
dni e n / e 22 

— = -Ki + 9Ki, — = (2 + 0'*2 - «1 + 16kiK2)K2 

^ = (| + 9k^)'«3, ^ = 2(1 - - 4«? + 4kiK2, (6) 

where e = 8 — rf. Above d = 8 the Gaussian fixed point is 
stable, while for d < 8 a runaway flow to strong coupling 
occurs, reminding of the one observed for the quantum 
phase transition of metallic models 1^ and of transverse 
field models |^ . The runaway flow is expected within the 
first order regime, u'--^"' < 0, but a strong coupling tricrit- 
ical fixed point limiting the known second order regime 
is still to be found. 

The RG for the DIC j/4 showed that its long-distance 
behaviour is dominated by a K^~contribution (like in | [Tl| ] 
but) for di"^ = 8 < d < 10. This leads to the modi- 
fied MF exponent 9^ = -j^, which satisfies the scaling 

relation 63 = 2/ /3a3 in d'"^^ = 8 and reduces to the MF- 
result in 10 dimensions. The dimensional shift by 2 in 
comparison with Ref. ||ll[ is due to coupling t. 

For given nonrandom chemical potential the frustrated 
spin interaction generates weakly nongaussian statisti- 
cal fluctuations of the fermion filling and vice versa, 
half filling exempted. Imposing instead a Gaussian 
distribution one finds the present problem mapped onto 
onto the Q-static approximation |]l3| of a metallic Ising 
spin glass in the limit of electron hopping range zero. 
This limit turns random hopping matrix elements into 
random site-local energies, which are equivalent to the 
fluctuating chemical potential and render the metallic 
Ising spin glass classical and static. The main new effect 
of /i-randomness is the generation of a classical T = 0- 
transition at < (Sn)^ >c= for < fi >= 0. 

We find that the random-/^ /S'G'/-model and the 
metallic Ising spin glass show tricritical behaviour on the 
T = 0-axis as well as in their thermal transitions with 
discontinuous low- and high-filling-regimes. At T = 
we find the tricritical point of the metallic spin glass with 
gaussian random hopping at 



where 2Eo denotes the width of the semielliptic elec- 
tronic band. Quantum dynamical corrections can be 
approximated by a generalized Miller-Huse method 0| . 
Details of the present work will be given elsewhere. 
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